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Abstract 

We introduce the peak normal form for elements of the Baumslag- 
Solitar groups BS(p, q). This normal form is very close to the length- 
lexicographical normal form, but more symmetric. Both normal forms 
are geodesic. This means the normal form of an element u~ l v yields 
the shortest path between u and v in the Cayley graph. For horocyclic 
elements the peak normal form and the length-lexicographical nor- 
mal form coincide. The main result of this paper is that we can 
compute the peak normal form in polynomial time if p divides q. 
As consequence we can compute geodesic lengths in this case. In 
particular, this gives a partial answer to Question 1 in [3]. 

For arbitrary p and q it is possible to compute the peak normal 
form (length-lexicographical normal form resp.) also the for elements 
in the horocyclic subgroup and, more generally, for elements which 
we call hills. This approach leads to a linear time reduction of the 
problem of computing geodesies to the problem of computing geodesies 
for Britton-reduced words where the i-sequence starts with t -1 and 
ends with t. 
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1 Introduction 



Baumslag-Solitar groups were introduced in [T] and they enjoy many 
remarkable properties, see e.g. Lyndon and Schupp [9\. The Baumslag- 
Solitar group BS(p, q) is a one-relator group defined by 

BS(p,g) := (a,t \ taPr 1 = a q ). 

The word problem is decidable in linear time, but it is still not 
known how to compute the geodesic length of elements efficiently 
Polynomial time algorithms for this problem were known only for hor- 
cyclic elements, [5j [6] or in the case where HS(p,q) is solvable, i.e., 
the case where p = 1, see [3]. As usual, a horocyclic element is an 
element of the subgroup (a) generated by a. More precisely, the paper 
of Murray Elder [3] presents a linear time algorithm how to compute 
geodesies for all words, when BS(p, q) is solvable. 

The result of Elder is the starting point for our paper and we 
generalize his result to the case where p divides q, thereby giving a 
partial answer to Question 1 in [3]. 

For this purpose we introduce the notion of peak normal form. 
The peak normal form is geodesic. Thus, it represents an element by 
a geodesic word in the Cayley graph of BS(p, q). There can be ex- 
ponentially many different geodesies, the peak normal form chooses a 
unique one; and it is defined by a natural condition. For horocyclic el- 
ements the peak normal form and length-lexicographical normal form 
coincide. 

Our main result is Theorem 17. H it states that peak normal forms 
can be computed in quadratic time, if p divides q. Actually, we have a 
more precise result, which yields Elder's linear time bound for p = 1. 

Our technique relies on the fact that length-lexicographical normal 
forms of horocyclic elements can be computed in linear time. See also 
[31 0(6] for similar approaches. We also extend the linear time result to 
elements which we call hills. These are words whose Britton reduction 
can be written in the form 

f3\t ■ ■ ■ Pktaot~ 1 ai ■ ■ ■ a(t~ x 

where all on and (3j are horocyclic. Moreover, we give a linear time 
reduction of the problem of computing geodesies to that of computing 
geodesies for Britton-reduced words where the t-sequence starts with 
t~ l and ends with t. 
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2 Notation and preliminaries 



Throughout the paper, let 1 < p < q be fixed positive integers. By 
a,di, /?, /9j,7,7j,5, (J,,v, p,cr,T we always mean integers. We reserve 
special fixed letters t, T, a, and A. By 6, Q{ we mean either t or T. 

The alphabet { t, T, a, A } is ordered by putting t < T < a < A, 
and it is equipped with an involution by a = A, A = a, t = T, and 
T = t. 

The involution is extended to words by 

a± ■ ■ ■ a m = ■ ■ ■ a\ for a, S { t, T, a, A } . 

We read A = a" 1 and T = t^ 1 in the Baumslag-Solitar group: 

BS(p,?) := (a,t | to^t" 1 = a q ). 

Every word w € {t,T, a, A}* can be read as w £ BS(p, q) and, of 
course, u; = w~ 1 6 BS(p, g). 

The Baumslag-Solitar group BS(p, g) is an HNN-extension of Z = 
(a) with stable letter t, and the HNN-extension is defined by the 
canonical homomorphism between the subgroups pTL and qL mapping 
the subgroup generator p to q. The group BS(p, q) can also be de- 
fined by a confluent and terminating string rewriting system BS over 
the alphabet { t, T, a, A } with the following rules, where denotes the 
empty word: 

aA 
tT 
a q t - 
aPT — 

Termination is not completely obvious, but the proof is standard with 
string rewriting techniques as explained e.g. in the textbook [7j. 
The rewriting system BS defines a congruence relation 

{ t, T, a, A }* x { t, T, a, A}* . 

BS 

This gives normal forms by computing w ==> w such that w is ir- 

BS 

reducible. Since BS is confluent, two words u and v are equal in 
BS(p, q) if and only if u <^=> v if and only if there is some w with 

BS 

u =^>- w <== v. Moreover, if w is irreducible, it is uniquely defined by 

BS BS 

U. 

In the following we write u ~ v if u <^4> v. Thus, u = v means 
identity as strings, whereas u ~ v means equality in BS(p, q). 



— ► 


Aa — 


-» 


— ► 


Tt — 


-> 




At - 


-» a q ~HA p 




AT - 


-» aP~ l TA q 
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The problem is that the length of w can be exponential in the 
length of w. To see this, consider e.g. p = 1 and the word t n aT n . Its 
irreducible descendant has length q n . 

The purpose of the system BS is therefore mainly to provide simple 
proofs for all basic properties about Baumslag-Solitar groups, see also 
[2] for more background about this approach. 

We investigate the problem of computing geodesies, i.e., given a 
word w over the generators a and t and their inverses A = a" 1 and 
T = t , find a shortest word g(w) over the alphabet { t, T, a, A } such 
that w and g{w) describe the same element of the group BS(p, q). 
The geodesic word g(w) is not unique, but its length |<?(ii>)| S N is 
well-defined. It is called the geodesic length of w. 

There might be exponentially many different geodesic words g{w) 
for a word w. So we are interested in unique normal forms as well. 
The length-lexicographical linear order is defined on words u and v 
by letting u <u v if either \u\ < \v\ or first \u\ = \v\ and second 
u is not behind v in the lexicographical order, which is defined by 
t <T < a < A. 

The length-lexicographical normal form of a word w is denoted 
by llnf(w). It is the first word v in this order satisfying v ~ w. 
Obviously, llnf(tf) is geodesic. Later we will introduce another (and 
more symmetric) geodesic normal form: the peak normal form. 

Confluence and termination of the rewriting system BS imply that 
the word problem for BS(p, q) is solvable. Therefore computing geo- 
desies (or any decidable geodesic normal form) is possible in finite 
(though using this naive approach at least exponential) time. On 
input w just enumerate all words up to length and check if they 
are equivalent to w. 

Throughout, we identify the word a a with the integer a, and we 
identify A a with —a. As a consequence, the free product Z*{ t, T }* of 
Z with the free monoid {t,T}* becomes both a quotient monoid (by 
a i— > +1 and Ah-1) and at the same time a subset of { t, T, a, A }*, 
where denotes the empty word. For elements of this subset Z * 
{ t, T }* we use normal forms. In other words u £ Z * { t, T }* is 
always represented as a sequence which alternates between integers 
and elements from { t, T } + . 

More precisely, we identify Z * { t, T }* with the set of sequences 
of words which alternate between words in { a }* U { A }* and letters 
in {t,T}. Thus, (42, t, t, 3, T, t) is allowed, but (43, — 1, t, t, 3,T, t) is 
not, although they denote the same element in Z * {t, T}*. (This 
means Z * { t, T }* becomes the regular subset of { t, T, a, A }*, which 
is defined by forbidding factors a A and Aa.) 

Note that mapping a word w £ { t, T, a, A }* to the corresponding 
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word u € Z * { t, T }* does not increase the length. It is also clear that 
all geodesic words in { t, T, a, A }* belong to the subset Z * { t, T }* C 
{t,T,a,A}*. 

Every word «£Z*{t,T}* can uniquely be written as a sequence 
u = a 9iai ■ ■ ■ 6 k ak 
with k > 0, ctj € Z, and B; L € { i, T }. Its length is 

k 

\u\ = k + y~] | Oi| . 

i=0 

Let || a || denote the geodesic length of an integer a. We define the 
norm \\u\\ of u by 

ll n ll = k + ll Q ill ■ 

4 = 

Let |<?(w)| be the geodesic length of u. Then we get 

k k 

\g(u) \ < \\u\\ = k + \\ a i\\ — \ u \ = k + ^2 l a *l ' 

i=0 i=0 

having equalities when u is geodesic. 

The objective is therefore to compute on input u a word v which 
minimizes ||u|| among all words v with u ~ v. In the next section 
we show that, given a word, we can compute an equivalent Britton- 
reduced word in polynomial time. Thus, in order to compute geodesies 
we only have to consider Britton-reduced words. 

Observe that ||a|| = |a| for an integer a implies |a| < 3q, and as 
soon as |a| > 2q, we find a geodesic using letters t (and T). This 
is trivial, because [iq ~ tfipT. It is also well-known (and shown in 
Section H|) that, if there exists some geodesic using the letter t, then 
there is some geodesic g\{a) which starts with the letter t, and there 
is some geodesic 52(a) which ends with the letter T. 

As a consequence, let w be a word and let u = ao^iai • • • @k a k 
minimize ||it|| for all u ~ w. Consider some < i < k where ^+1 = t. 
No geodesic of can end in a T, so |<?(aj)| = \ai\ < 2q. As we have 
(±g)i ~ t(±p) we see that actually |aj| < q. The same happens if 
1 < j < k and 8j = T, then |<x,-| < q, too. 

Thus, large values |aj| can be found at local peaks of u, only. A 
local peak of u is a position i G { 0, . . . , k } such that < i < k and 
9i T and #,,+1 7^ t. Note that is a local peak for k = 0. For k > 
it is a local peak if 9\ = T. 

The notion originates from the shape of the path of n, if reading 
a T means going one step downwards and reading a t means going 
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one step upwards, as illustrated by the following figure for the word 
u = ATa 2 taTattA 10 TaTa 2 tta A2 T: 





-10 


42 


-1 


1 / Yl 




\ 


/ v \ 


I 


A word with four local peaks 



As integers aj with large absolute value \o.i\ can only be found at 
local peaks of u, the number of local peaks is an important parameter. 
However, for technical reasons which will not become apparent until 
later, we prefer to count the number of sinks. A sink is dual to a local 
peak. Formally, a sink of a word u = ao^i«i • • • &k a k is a position 
i E { 0, . . . , k } such that 9i ^ t and ^ T. Note that is a sink 
for k = 0. For k > it is a sink if Q\ = t. The difference between 
the number of local peaks and the number of sinks is bounded by ±1. 
For example, if p = 1 (which means that BS(j>, q) is solvable) Britton- 
reduced words may have 1 or 2 local peaks, but always exactly 1 sink. 
Here are the pictures for words having one sink and one or two local 
peaks: 
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a 2 ta 13 



-1 42 



2+„42 




ATaHa 



We emphasize this is in a special notation: Let u = ao^iai • • • #fcOjfc 
be a word. The number of sinks is denoted by s(u) which is defined 

as 

s(u) = \{i G { 0, . . . , k } I i is a sink }| . (1) 



3 Britton reductions 

Let u = ao#iai • • • Qk^k be a word in Z * { t, T }*. A Britton reduction 
step means replacing in u = ao#iai • • • Ok®k either a factor tfipT by 
/iq or a factor T^iqt by /ip, where /i G Z. After that, we rewrite the 
new sequence by the corresponding word without factors aA and Aa 
in the free product Z*{t,T}*. So a Britton reduction step decreases 
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the number of letters from { t, T } but it may increase the length. If 
e.g. 9iOi6i + i = tfxpT, then one Britton reduction step yields 

ui = a 9iai ■ ■ ■ 9i-i(ai-i + fiq + Oj+i)6>i+ 2 ai+2 • • • 9 k a k . 

A Britton reduction may increase the length by a factor q/p, but 
the important point is that < Indeed the geodesic length 

of Oj_i + fiq + Oj+i is at most |aj_it^pToj+i| . A Britton reduction 
neither increases the number of local peaks nor the number of sinks. 

A word is Britton-reduced, if no Britton reduction step is possible. 
As the coefficients Oj can increase exponentially, from now on we keep 
all integers in binary notation. The number of bits remains linear in 
the input length because the sum over log 2 (|aj|) is never greater than 
q\w\. For complexity considerations it is convenient to assume that 
an arithmetic operation on n-bit integers (with n € can be 

performed in constant time. Then each step in a Britton reduction 
needs constant time. So we can produce a Britton-reduced equivalent 
word (with integers written in binary) in linear time. The number of 
bit operations is actually not worse than quadratic (if polylog terms 
are ignored). 

Britton-reduced forms are far from being unique. On the contrary, 
there can be exponentially many Britton-reduced equivalent words. 
But if w ~ u = oo#iOi • • • ^fcOfc and u is Britton-reduced, then the 
t-sequence 

(0i,...0 fc ) e {t,T}* 

is uniquely determined by the word w. (This can easily be seen with 
the help of the rewriting system BS defined above since applying its 
rules leaves a Britton-reduced word Britton-reduced.) 

Thus, as a Britton reduction never increases the norm, we can 
assume that the input word u = Qo^iQi • • • 9 k a k is Britton reduced and 
the integers Oj are written in binary. The objective is therefore reduced 
to minimizing the norm ||n|| (and to computing geodesic normal forms) 
for Britton-reduced words and to computing geodesies of horocyclic 
elements. 

As usual, a word u £ { t, T, a, A }* is called horocyclic, if its im- 
age in BS(p, q) belongs to the horocyclic subgroup (a) generated 
by a. We adopt this notion here. Thus, a Britton-reduced word 
u = oo^ioi • • • #fcOfc is horocyclic if and only if k = 0. As a con- 
sequence, membership of the horocyclic subgroup can be tested in 
linear time. 

The next step is to compute geodesies of horocyclic words. 
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4 Geodesies for horocyclic words 



Techniques for computing geodesies for horocyclic words can also be 
found elsewhere, see [3j[5j[6]. Corollary |4.7l and Remark 14.81 have been 
inspired by recent results by Freden et al., see [HI Prop. 9.1]. Since 
our presentation is however quite different from other approaches we 
give full proofs. 

The strategy to compute length-lexicographical normal forms for 
horocyclic words w ~ a G Z relies on the following simple observa- 
tions: 

1. ) Horocyclic elements commute in BS(p, q). 

2. ) If \a\ is small, then we can use table lookup. 

3. ) If \a\ is large, then llnf(tt;) begins with t, see Lemma 14. II 

4. ) If w ~ tuTv, where the indicated letters t and T match in a 

Britton reduction, then w~(3€Z and a = (3 mod q. 

5. ) If w = uw'v and llnf(u>) = uw"v, then llnf(u/) = w" . 

We begin with a lemma which might be of independent interest. 

Lemma 4.1 Let w ~ a be horocyclic and g(w) be a geodesic repre- 
sentation. Then we can write: 

g(w) = (3 k t--- (3ita Ta± ■ ■ ■ Ta k . 

Moreover, there are also geodesic representations gi(w) and g2{w) 
such that 

gi(w) = t- ■ -ta r(ai + • • • T(a k + f3 k ) 
g 2 (w) = (a k + k )t • • • (ai + f3i)ta T ■ ■ • T 

We also find such geodesic representations with k > 1 as soon as 
\a\ > 2q. 

Proof. The shape a ~ g(w) = (3 k 0' k ■ ■ ■ /3i^ao^iai ■ • • 9 k a k is 
obvious because in the Britton reduction all t and T must vanish. 
Clearly, a 2q ~ ta? p T and A 2q ~ tA 2p T, so the assertions are trivial for 
k = 0. Let fe > 1. 

1.) Let g(w) = uv be the product of shorter horocyclic words. If 
both u and v have geodesic representations using the letter t 
(or T), then we can write g(w) ~ g2{u)g\{v) with | | = 
\92{u)gi{v)\- But g2(u)gi(v) is not geodesic, because the first 
letter of g\{v) is t and the last letter of g2(u) is T. Since 
g(w) ~ vu ~ uv, we may assume that v = a k . The claims 
follow by induction, because elements of Z, such as a k and (5 k 
commute with horocyclic elements. 
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2.) If g{w) is not the product of shorter horocyclic words, then we 
have Pk = ak = and either g(w) = Tut or g(w) = tuT. If 
g(w) = tuT, then we are done by induction. So assume by 
contradiction g(w) = Tut. 

i. ) Let k > 2 and g(w) = Tut. This is impossible, because 

k > 2, therefore g(w) ~ Tgi(u)t, and gi(u) has t as its first 
letter. 

ii. ) Let k = 1 and g(w) = Tat. Then we have a G gZ, so = 

Tfiqt for some G Z. But T/igt ~ /xp and \/j,p\ < \/j,q\ + 2. 
Thus we find a contradiction again. 

□ 

Example 4.2 Lei p = 1 cm<i q = 2. Our techniques will show that 
t k a(Ta) k is a geodesic representation. The word is obviously horo- 
cyclic. It is also clear that Lemma \^.l\ can be generalized such that in 
this case we find exactly 2 k geodesic representations: 

k k 

([[a £ H) a (fjTa 1 " 5 ') where e { G {0, 1}. 

i=l i=l 

The first phase in our computation to produce the length-lexi- 
cographical normal form for horocyclic elements will a greedy linear 
time reduction to compute the length-lexicographical normal form of 
Britton-reduced words where all 9i are equal to T. For this purpose 
we introduce the following notion. A slope is a word w which has the 
form 

w = aoTai ■ ■ ■ Tak- 
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The slope a 42 TA 3 Ta 5 TTA 



In [3] slopes are called words of type N. 
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Proposition 4.3 Let w £ {t,T,a, A}* be a horocyclic word with w ~ 
a£Z, Then we can compute in linear time a slope u and a natural 
number I £ 0(log \a\) such that 

llnf (w) = t e llnf (u) 

Proof. We begin with some precomputational steps. 

1. ) Replace w by the corresponding word u = <5o#i<5i • • • 6 m 5 m in the 

free product 7L * { t, T }*. 

2. ) Compute a Britton-reduced equivalent word u' = S'qO'^ ' " " Q't'&l' 

keeping integers in binary representation. If £' ^ 0, then w was 
not horocyclic and we can stop. 

3. ) Now we may assume that w = a £ Z. 

Next, we perform a greedy linear time algorithm. By symmetry 
we assume < a £ N. The basic idea is that as long as a > 2q, the 
unary notation a a is not a length-lexicographically first representation, 
because ta 2p T <u a 2q . This motivates rewriting a = q/j, + with 
< < q. We can replace a by tpfiTfi, and we have < p[i < a. 
Repeating this greedy process for \xp as long as pfi > 2q yields in 
linear time a representation w ~ t l 0oT0\ ■ ■ ■ T0g with £ £ 0(log |a|), 

< 0i < q for i / 0, and < O < 2q. 

Define Oj = t % 0§T0\ ■ ■ ■ T0%. For each i we found /Xj > such that 
cti = qui + 0i and 

fpoT-'-pi-tT ~ Pfli . 

We have a = ae and is enough to show that the word llnf (a') begins 
with t % for each a' with a, < a'. This is trivial for i = 0. Now let 

1 > 1. We know 2q < cti < a', so by Lemma |4. II we obtain: 

llnf (a') = f'joTPi ■ ■ ■ TP\, 

for some i! > 1. Moreover, /?^ ; < q. We have a' = 0[t mod g. Hence 
there exists fi' > 1 with a' = //g + 0' v Recall that oti = ^q + /3j where 
< Pi < q- Therefore on < a' implies ^ < fi'. Induction applies, and 
we can conclude that llnf(p^') begins with t l . Thus, llnf (a') begins 
with t\ □ 

By Proposition 14.31 it is enough to compute length-lexicographical 
normal forms for slopes. Essentially, this can be performed by a finite 
transducer. Our algorithm relies on dynamic programming and needs 
some preparation. 

Let R denote the set (of constant size) defined by 

R = {p£Z | \p\ <r }, 
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where in this section the constant r is the least positive integer such 
that 

r + q — 1 

r>p 2 +g-l. 2 

Actually, minimizing the concrete value of r is of little importance. 
All we need is that for all constants c > there is a positive integer 
r = r(c) such that 

r + c 

r > p ■ h c. (3) 

So, in Eq.[2]the constant c was chosen to be q—1, which guarantees 
in particular that < q implies (3 € R. 

Proposition 4.4 There is a linear time algorithm (where the number 
of bit- operations is linear, too) which on input a slope w = fioTfii ■ ■ ■ T@£ 
with 

|#| <q fori^O, and \/3 \ < 2q (4) 

computes the length-lexicographical normal form llnf(tt;). Moreover, 
we have 

llnf (w) = t k (3T ai ■ ■ ■ Ta k+e 

with k G 0(1). 

Note that the slope u given by Proposition 14.31 satisfies (j3J). Thus: 

Corollary 4.5 The length-lexicographical normal form of horocyclic 
elements can be computed in linear time (where the cost of arithmetic 
operations is considered to be constant). 

Proof. (Proposition I4.4j ) The input w = (3qT(3\ ■ ■ ■ Tfin is Britton- 
r educed. 

We initialize a table with entries llnf(/o) for all p G Z with \p\ < r+q 
in constant time. The desired form of the entries is guaranteed by 
Lemma 14.11 For i = we find the information in this table because 
|/?o| < 2q. Now, let t > 0. 

For each < i < i and each 7 £ R we define the word u(i, 7) by 

u(i,j) = PoT-'-Pi-iTj. 

We show that we can compute llnf (u(i, 7)) efficiently by table lookup. 
Recall that f3 e € R. 

If i = then 14(0,7) = 7 is a small integer, and Unf(7) is in the 
precomputed table. 

Let i > and assume for simplicity (and for a moment) that for all 
7 £ R the values llnf (u(i, 7)) have been computed and stored in a table 
of size \R\. We wish to compute the length- lexicographical normal 
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form of u(i + 1, pj+i) = [3qT ■ ■ ■ j3iTpi + \ for each pj + i G i?. The word 
llnf (u(i+l, Pi+i)) ends in some Tj, where I7I < q and pi + \ = 7 mod g, 
because a slope is Britton-reduced. Thus, pi + \ = pq+j for some p € Z 
and we obtain Tpi + \ ~ ppTj. Now we have |pj+i| < r. Hence for 
i = we obtain 

r + g — 1 

MP + Pi < P ■ \-2q - 1 < r + q. 

q 

For i > we obtain 

r + g — 1 

pp + « < p • h Q — 1 < r. 

Thus, in both cases we can compute by table lookup: 
llnf(-u(i + 1,/Oj+i)) = min{llnf(n(z,pi))r7} , 
where the minimum is taken over all pi £ R and | —y | < q such that 
3p : pi = f3i + pp A pi + i = pq + 7. 

So far we have proven Proposition 14.41 except for the linear time 
bound. In fact, we have shown that the lengthdexicographical normal 
form of slopes (and thus of horocyclic elements) can be computed in 
polynomial time. It remains to explain how we manage to find the 
minimum in constant time with a constant amount of information. 
Unfortunately, the explanation is slightly technical. If the reader is 
interested in polynomial time results only, he or she is invited to skip 
the rest of the proof. An example of how the method works can be 
found in the appendix. 

For the explanation we observe first that 

I |llnf(u(i,p))| - |llnf(n(i,r))| | < 2r 

for all p, t € R. One idea is therefore that instead of keeping the 
words Unf(n(i,p)) in the table, we store in it only certain suffixes of 
these words of maximal length 2r. 

This would be even more evident if it were enough to compute the 
geodesic length. Then the table would just need to store the length 
differences 

\g(u(i,p))\ - \g(u(i,0))\ e {-2r,...,2r} 

in addition to the absolute geodesic length of one of these elements, 
say, \g(u(i,0))\. 
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As we are more ambitious than simply computing geodesic lengths, 
we need more subtle data structures. More precisely, factorize each 
llnf (u(i, p)) as 

llnf(u(i,p)) =Pi, p • s i>p 

such that \pi )P \ = min{|llnf (u(i, r))| : r £ i?}, thus having \si tP \ < 2r 
for all p £ i?. We modify the above algorithm so that after the i-th 
round we only the store the following information: 

• A list of the Si >p for p £ i?, and 

• the lexicographical ordering of the prefixes pi iP (p £ R). 

Note that this information only takes a constant amount of space and 
that all pi tP have the same length. 

It is clear now that we can perform the minimum search in constant 
time. Say, we need to compare Hnf(u(i, p))T*y and llnf {u(i,r))TS. As 
\Pi o\ = \Pir\ we can compare the length by looking at the lengths of 
Si^pT^j and Si^ T T5. In case \s^ p T~(\ = \si tT TS\ we check the ordering 
between pi tP and pi jT . If p^ p = pj jT , we finally compare Si :P T^ and 
Si, T T5. 

The length of the shortest word llnf (it (z, p)) (with p € R) strictly 
increases from round i to round % + 1, since a factor of type T7 is 
concatenated. So we can actually update the information in constant 
time. We simply have to compute new suffixes and a new linear order 
on R. As we need to recover Pi+i tP later, the algorithm outputs a 
column vector with the portions of the suffixes which were cut due to 
length constraints plus a pointer to the index of the preceeding column 
where the minimum was achieved. 

The output of the algorithm can be viewed as a matrix with 2r + 1 
rows and a linear number of columns. The entry (p, i) contains a word 
of constant length and a pointer to some (r, i — 1). In the last round I 
the output consists of the suffixes S£ tP and, using these, the length-lexi- 
cographical normal form for each u(£, p) can be read by a single scan 
from right to left in the matrix, starting at position (p, t) and following 
the pointers. In particular, we can read off llnf(u>) = llnf(n(£, (5^)) in 
linear time. □ 

Corollary 4.6 Let w £ { t, T, a, A }* be a horocyclic word with w ~ 
a £ Z and let I £ N be the number computed in Proposition ^. 3\ Then 
we have 

llnf (w) = t k (3T ai ■■■Ta k 

with < k — I £ O(l), |ci!j| < q, = a mod q, |/?| < 2q, and if k > 
then (3 £ pZ. 



13 



Proof. Propositions 14.3 1 and 14.41 show that llnf(u)) has the form 
t k [3Ta\ ■ ■ ■ Tctk with k — I € 0(1). The other assertions are a direct 
consequence: \ai\ < q, = a mod q, \{3\ < 2q, and if k > 0, then 
(3 € pL. □ 

Let us call a slope w = PoTfii ■ ■ ■ Tfy semi-horocyclic, if there is a 
number k such that t k u is horocyclic. Clearly, if such a /c exists, then 
we have k = I, so k is unique. 

Corollary 4.7 

1. J Tie set of slopes in length-lexicographical normal form 

{ w = faTfr • • • Tf3g G {T,a,A}* | w = llnf (w) } 
is regular. 

2. ) If p devides q, then the set of semi-horocyclic slopes in length- 

lexicographical normal form 

{ w € { T, a, A }* | w = llnf (w) and w is semi-horocyclic } 
is regular. 

3. ) If p devides q, then the set of horocyclic elements in length-lexi- 

cographical normal form 

{ w € { T, a, A }* | w = llnf(u>) and w ~ a G Z } 

is a deterministic (and unambiguous) linear context-free (one- 
counter) language; and it can be recognized in log-space. The 
growth series of the horocyclic subgroup is a rational function. 

Proof. Formally, Statement [T] does not follow from assertion 
in Proposition 14.41 but analyzing its proof shows that all arithmetic 
computations concern only a constant number of integers of constant 
size. This can be done in the finite control. For accepting a length- 
lexicographical normal form, we do not need any output, we just have 
to check that the input agrees locally with a potential output. Again, 
this can be done in the finite control. The result follows. 

Statement [2] follows from [T] because a test whether a slope is semi- 
horocyclic can be done by counting modulo p. 

Statement follows from [21 we simply have to check additionally 
that the number of t and T match. This can be done by a deter- 
ministic (one-counter) pushdown automaton, which translates to an 
unambiguous linear context-free grammar. The growth series of un- 
ambiguous linear context-free languages is rational, [8]. □ 
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Remark 4.8 Statement is essentially a result due to Freden et 
al, see J?J Prop. 9.1]. Statement^ is slightly more general than J(5j 
Prop. 9.1] because our result holds for all p and q. There is a crucial 
difference for p \ q. It is only when p \ q that we can test with a 
push-down automaton whether an input word w = t e (3oTf3i ■ ■ ■ Tftg is 
horocyclic, see J<3 Thm. 7.2]. So we need this as a promise in order 
to produce horocyclic words with the sort of finite transducer we used 
in the proof of Proposition \4-4\ A transducer cannot compare m = I; 



and even if m = i, it cannot test whether the input is horocyclic. This 
part is however trivial for p \ q by counting modulo p. This is why 
follows from [7] easily in this case. 



5 Peak normal forms 

We consider A = Z U { t, T } as an infinite alphabet of symbols with 
the following linear order: 

1. ) t < T < a 

2. ) a<p,i£ \a\ < \p\ 

3. ) a < 0, if < a = \(3\ 

The length-lexicographical order on A + transfers to a linear order <a 
on words of Z * { t, T }* C {a, A, t, T}* . 

Consider a word w € Z * { t, T }* such that 

W = UqU^i ■ ■ ■ ttpCtg. 

For every position % G { 0, . . . , I } of w we define its height by 
h(0) = and h(i) = h(i - 1) + 1 if B\ = t or h(i) = h{i - 1) - 1 when 
8[ = T. The height of w is defined by h(w) = max { h(i) | < i < t }. 

Again, this notion arises from the shape of the path of u, if reading 
a T means going downwards and reading a t upwards, as we have done 
before. Consider e.g. w = ATaHaTattA 10 TaTaHta 42 T: 



1 


-10 


42 


-1 

< 


i /V 




V 

1 \ 


/ V/ \2 




A word of height 1 



For the peak normal form we compute some Britton reduction first. 
Thus, we let w ~ u and 

w ~ u = ao^iQ^i ' • • QkOtk 
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where u is Britton-reduced. Recall that the sequence (6±, ... ,6k) de- 
pends on w only. 

We say that the position i is the peak of u, if i is maximal among 
all i with h(i) = h(u), i.e. it is the rightmost among the hightest local 
peaks. Define u\ by 

u\ = a 0i ■ ■ ■ aj_i0j_i 
We obtain a natural factorization 

U = UlOiU2- 



peak 




The peak of a word 



We say that a Britton-reduced word u = ao9±ai ■ ■ ■ Q^ctk is the 
Britton peak normal form of w if 

1. ) w ~ u and = 

2. ) Among all choices satisfying w ~ it and = ||tt|| we choose 

the one where u\ is the first in the order <a, after that we 
minimize «2 in the order <a- 

We say that a word v is the peak normal form of w if 

v = llnf (a o )0i llnf (ai) • • • 9 k llnf (a fc ) 

where u = ao$i a i ■ ■ ■ ^fcOfc is in Britton peak normal form. The peak 
normal form of a word w is geodesic, and it is denoted as pnf(io) in 
what follows. 

6 Difficult cases and solution for hills 

Let u = ao^iai • • • Ok a k be Britton-reduced. Due to Section H] we can 
compute the norm ||it|| and length-lexicographical normal forms for 
each ctj. Unfortunately, we are still not able to compute a geodesic 
for u efficiently, in general. But at least we can identify difficult cases 
and solve the the problem for so-called hills. A hill is a word w such 
that we have 

w ~ u = Pit ■ ■ ■ j3\taoTa\ ■ ■ ■ Ta m 
for some u and £, m > 0. 
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A hill 



Note that all horocyclic words are hills, so their length-lexico- 
graphical normal form is a hill representation. It is clear that w 
is a hill if and only if its Britton reduction already has the form 
fyt ■ ■ ■ PitaoTai ■ ■ ■ Ta m . 

We show that we can compute the peak normal form for hills w 
very efficiently, see also [3]. If w is not a hill, then we can in linear 
time reduce the computation of the peak normal form of w to the 
the computation of the peak normal form of so-called difficult words, 
which are defined below. 

For p | q we will solve the remaining difficult cases in the next 
section. 

In this paper we call a word w difficult, if its Britton reduction is 
ao#i«i ' • • QkOtk with 9\ = T and 9k = t, in other words it has the form 
aTvtf3. 




A difficult word 



Note that every Britton-reduced word u has a unique representa- 
tion 

u = a\t ■ ■ ■ attDTflm ■ ■ ■ Tj3\ 

where either D = 5 is horocyclic (making u a hill) or D = aTvtfi, i.e. 
D is difficult. 
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Factorization of a word 



Assume that we have an algorithm that computes peak normal 
forms for difficult words. Let Time(d) denote the maximal time to 
compute pnf(D) for ||D|| = d for = d. We assume that Time(d + 
c) € 0(Time(d)) if c is a constant. We will show later that this 
assumption is justified in the case p \ q. It remains reasonable in 
general, since it holds for every sensible complexity bound, such as 
polynomials or singly exponential functions. 

The following result generalizes Proposition 14.41 For its proof we 
use the constant r as defined in Eq. [2j 

Theorem 6.1 Let u = a^t ■ ■ ■ a\tDT 'f3\ ■ ■ ■ T '(3 m be Britton-reduced 
with D either horocyclic (D S 7L) or D difficult. 

1. ) If D is horocyclic, then we can compute the peak normal form 

pnf(u) in linear time. 

2. ) If D is difficult, then we can compute the peak normal form 

pnf(tt) in 

0(Time(||D|| + max{||a;|| , ||^-||}) + \\u\\). 

Proof. Applying a greedy algorithm similar to the one used in 
the proof of Proposition 14.31 we can ensure that < |o!j| , \ j3j\ < q for 
all 1 < i < k, 1 < j < m. Doing this, D might become longer, but 
this doesn't affect the case 1.) since we have a linear time algorithm 
for horocyclic words. In case 2.) it suffices to proove a time bound of 
0(Time(||Z>||) + |M|). 

Now let 

u(i,j, p, 5) = ptai^it ■ ■ ■ a\tDT(3\ ■ ■ ■ Tf3j-\T5. 

By induction on i + j we prove that the theorem holds for every word 
u(i, j, p, 5) with p, 5 £ R = {~{ e Z \ | -y j < r }. 
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If i + j = 0, we need to compute the peak normal form of pD5, a 
word of length at most \D\ + 2r. If D is horocyclic, Corollary 14.51 tells 
us that this can be done in linear time. If D is difficult, then we need 
0(Time(||D|| + 2r)) which, by our assumption on Time, is bounded 
by C(Time(p|[)). 

Now consider % + j > 0. By symmetry we may assume that i > 0. 
Since D and u(i,j,p,5) are Br itton-r educed, pnf (u(i, j, p, 5)) starts 
with jt and we have 7 = p mod q. In addition, (7! < q, since |a 9 i| = 
\A 9 t\ = q + l>p+l = \ta p \ = \tA p \. So, since the peak is inside of 
D, like in the horocyclic case, we have 

pnf(u(i,j,p,5)) = min{7tpnf(n(z - l,j,p',5)) } 

where the minimum is taken over all 7 and p' such that I7I < q and 
p = 7 + pq and p' = pp + a.i-\. Again, since p G R, we have € i?, 
so induction applies. 

For the implementation we use again a table of size \R\. It stores 
the length information 

\pnf(u(i,j,p,5))\ - |pnf(u(i,j,0,<5))| G { — 2r, . . . , 2r } 

as well as the ordering between various words pnf (u(i, j, p, 5)) and 
pnf(u(i, j,T, 5)) as a linear order on R. 

This allows the minimum search in constant time by table lookup in 
a table of constant size with constant size entries. Indeed the ordering 
between *yt pnf (u(i,j, p, 5)) and 5t pnf j, r, 5)) is dominated by the 
length and next by the ordering between -ft and 5t. For -ft = 5t we can 
refer to the ordering between pnf (u(i, j, p, 5)) and pnf(n(«, j, r, 5)). 

The table update is possible in constant time, too. The output is 
produced from right to left in this phase. □ 

Corollary 6.2 Let w be a hill. Then we can compute the peak normal 
form pnf (id) in linear time. 

The proof of the following Corollary 16.31 is rather technical. We 
do not use the result anywhere, so we leave its proof to the interested 
reader. 

Corollary 6.3 Let p divide q. Then the set of hills in peak normal 
form 

{w \ w = pnf(w) and w is a hill } 

is a deterministic (and unambiguous) context-free language. Rs growth 
series is a rational function. 



19 



7 Complete solution when p divides q 



Prom now on we assume that p divides q. We give an algorithm 
for computing geodesies of any element in BS(p, q) which runs in 
quadratic time, assuming that arithmetic operations in Z take con- 
stant time. One reason why things simplify is that we have a canonical 
mapping 

7T : BS(p, q) -> (Z/pZ) * { t, T }* 

which is induced by mapping a Britton-reduced word ao#iai • • • Qk a k 
to: 

(qo mod p)Q\{a\ mod p) • ■ ■ 9k{o<-k mod p). 

Using the confluent string rewriting system BS, we see that, in par- 
ticular, the sequence (cto mod p, . . . , mod p) € (Z/pZ) fc depends on 
the image of u in BS(p, <?), only. 

Recall that according to Eq. [1] the number of sinks is denoted by 
s(w). We have 1 < s(w) < 1 + ||w||. Theorem 17.11 is the main result 
of the paper. 

Theorem 7.1 Let p be a divisor of q. Let w 6 {t,T, a, A}* be an 
input word. Then we can compute a geodesic and its geodesic length 
in quadratic time 0(s(w)(s(w) + ||^||)). 

The rest of the paper is devoted to proving this result. The formal 
proof of Theorem 17.11 is postponed to Section 17.41 

We may start with a word w£Z*{t,T}* such that 

w ~ u = ao6\ai ■ ■ ■ QkOtk- 

Recall the definition of the height from Section [5] The word w is 
called a valley, if h(u) = and h(k) = 0. 




A valley 



A Britton reduction cannot increase the height and leaves the 
height of the last position invariant. Therefore, whether or not a 
word is a valley can be checked on its Britton reduction. 
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Valleys are generated by the following context-free grammar where 
S is an axiom: 

S — > SS | aTS(3t | a 

The following lemma can be based on this grammar. The proof is 
straightforward, but it uses in a crucial way that p \ q. 

Lemma 7.2 Let p \ q. Let v be a valley. Then pv ~ vp. 

Proof. We have p + a = a + p in Z. If u = vw and v, w are 
valleys, then, by induction, pvw ~ vpw ~ vwp. If v = aTw[3t, then 

pv = (p + a)Tw(3t 
~ aTqwf3t 

~ aTw((3 + by induction, because p | g 

~ aTwfltp. 

□ 

7.1 Reduction to valleys 

Remember that we assume p | g and it only remains to compute 
peak normal forms for difficult Britton-reduced words u = aTvt(5. 
Let I = h(u). There is a unique m € N such that T e ut m is a val- 
ley. The word T l ut m is still Britton-reduced! By Corollary 17.61 we 
can compute the peak normal form pnf(T^ni m ) in quadratic time. 
The word pnf (T^ut" 1 ) has the form a\T ■ ■ ■ a{Fw. Now, the sequence 
ir(T e u) begins with T l G (Z/pZ) * {t,T}*. In other words, we obtain 
a\ = ■ ■ ■ = ct£ = mod p. Each of the first I positions finds a match- 
ing position greater than I of the same height. So, we can shift all 
integers ctj for 1 < i < I to the right by Lemma [7.21 without increasing 
the length. Hence, the peak normal form begins actually with T . For 
m = we conclude that pnf (T e ut m ) = T e pnf (ut m ) = T e pnf (u) and 
we are done. 

It remains to compute pnf (it) for m > 0. In this case we write 

pnf(T^ut m ) = U17U2 = T l v\^ui 

where 7 is the integer on the (rightmost) peak of T e u, which is also 
the (rightmost) peak of u. The peak of u is not the last position in u, 
because m > 0. Hence, U2 is a valley which begins with the letter T. 
Moreover, we may assume that either u\ = or it ends with t. 

Now we concentrate on the valley 7^2 ■ This time we compute, 
again in quadratic time the peak normal form from right to left. This 
means we compute 

pnf(7U2") = 5v 2 t m , 
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where V2 begins with T. 

We claim that pnf(u) = v\5v2- Indeed, let pnf(u) = v^S'v^ where 
v 2 begins with T, ir(v2) = n(v 2 ) and either v\ = or v\ ends in t. 

Obviously, Vi5v2 ~ v'-yb'v' 2 and l^i^l = l^i^'?^!- Now, v[ <a v\, 
hence T^v[ <a T^v\ and therefore v[ =v\. 

As we have |i>i(k>2| = ki^'^l an< ^ v i = Vl we mus t have v 2 <a V2 



and v' 2 t m <a V2t m . Now, 5v2 ~ 5'v' 2 and hence 
We conclude 5v2 = 5'v' 2 , and finally 5 = 5'. 



7u 2 



5'v' 2 f 




7.2 Standard valleys 

A word V is called a standard valley, if V = ol\0\ ■ ■ ■ ockQkfl is a Britton- 
reduced valley such that (3 = 0, |ct!j| < q for all 1 < i < k, and < p 
whenever 1 < i < A; and 0i = T. 

Note that non-trivial standard valleys end in the letter t. By Sec- 
tion 17.11 we may assume that we start the computation of a peak 
normal form with a Britton-reduced valley. The next lemma is a key 
step for the sequel. 
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Lemma 7.3 Let p \ q and v be Britton-reduced valley, then we find 
in linear time a standard valley V and an integer 7 € Z such that 
v ~ Vj. 

Proof. First note that is always Britton-reduced, because a 
non-trival standard valley ends in the letter t. 

1. ) For v = a £ Z we can choose V = 0. 

2. ) Let v = uw where u, w are shorter valleys. By induction, there 

are a standard valley U and (5 € Z such that tt ~ and is 
Britton-reduced. Hence v ~ U (3w and is a Britton-reduced 
valley. By induction, there are a standard valley W and 7 G Z 
such that v ~ U (3w ~ UW^y. 

3. ) Let u = aTu(3t. Then u is a Britton-reduced valley; by induc- 

tion, there are a standard valley £/ and 5 € Z such that u ~ US. 
Write a = fip + a' with < a' < p, then we have 

v ~ q'T[/((5 + w + 

Write 5 + fj,q + (3 = vq + (3' with < /?' < q, then we have 

v ~ n'TUfltvp. 
The word = a'TU (3't is a standard valley. 

□ 

By Lemma 17.31 it is enough to compute geodesic normal forms for 
words u = Vw where V is a standard valley and w ~ 7 is horocyclic. 
Again, this can be done with a dynamic programming approach as we 
demonstrate in the next subsection. 

7.3 Computing geodesies for Vw 

In this section V is always a standard valley and w is word such that 
tw is Britton-reduced. 

Every standard valley can be generated by the following context- 
free grammar: 

S -> SS I aT5/3t I 0. 

Here S is the axiom and a, (3 denote intergers which are also viewed 
as terminal symbols and satisfy \a\ < p and \/3\ < q. Actually, the 
grammar produces non-standard valley as well, because we can pro- 
duce valleys which are not Britton-reduced. But this is not important 
for what follows. 



23 



The number of sinks (c.f. Eq. [T]) of a standard valley V admits a 
nice recursion: 



1 if V = 

s(V) = I s(U) + s(W) if V = UW and U + + W 
s{U) if V = aTUt/3. 

We now choose a constant r similar as in Eq.[3J but this time using 
the constant c = 3q — 2. More precisely, we let r G N be minimal such 
that 

r + 3g - 2 

r>p- . (5) 

The next step is to define for each standard valley V a range 
R(V) C pZ such that |JR(V)| G 0(s(V)). The precise definition and 
necessary properties of R(V) are determined by the next rather tech- 
nical lemma. 



Lemma 7.4 Let p \ q and let r be the constant of Eq. [H Given a 
standard valley V as input, we can (define and) compute in linear 
time a range R(V) such that the following two properties hold: 

1. ) R(V) Q{p£pZ \ \p\<r- s(V) }. 

2. ) For all p G R(V) there is a standard valley V p with V ~ V p p. 
Proof. 

1. ) For V = define R(0) = { }. 

2. ) Let V = UW and U ^ ^ W. We define 

R(V) = {p£Z | p = a + r,a £ R(U),t G R(W) } C pZ. 

Clearly, p G i?(V) implies |p| < r • (s(C7) + s(W)) = r • s(V). 
Let a G i2(C) and r G R(W). As a € pTL we can write 

V ~ U(jO~W ~ E/^ T (<T + r). 

3. ) Let V = aTU fit. For a G i2(C7) consider in a first step all values 

a 1 with a = ep + a' and \a'\ < p. We know [e| < 1. We can write 

V ~ a'TU a {a + eq + P)t. 

Note that |<t + eg + /3| < r • s(V) + 2g — 1. Thus a leads to some 
values (3' and hence some values p such that er + eq + (3 = p,q + (3' 
and j < q and p = pp. Note that 

r • s{V) + 3g - 2 
|p| = |p| P < P ■ < r ■ s{V ). 
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We define R(V) to be the set of all p which are possible outcomes 
for some a G R(V). 
Of course, we can write 

V ~ a'TU a (a + eq + 0)t ~ a'TU a [3'tp. 

□ 

Theorem 7.5 Lei p j g. VKe can design a quadratic time algorithm 
running in time C?(s(V)(|| V|| +s(V))) which solves the following prob- 
lem. 

Input: A standard valley V 

Problem: Compute for all p G R(V) words V p in peak normal 
form such that the following condition is satisfied. 
For all w, where tw is Britton-reduced, it holds: 

pnf (Vw) = min { V p pnf (jnv) | p G R(V) } (6) 

Proof. Note that for standard valleys the length-lexicographical 
ordering coincides with ordering in peak normal form. Thus, if we wish 
to test later whether pnf(C/ cr ) < pnf(U p ), it is enough to remember the 
length-lexicographical ordering between prefixes of the same length, 
and in case they are equal we can compare suffixes. During the proof 
we will point out where this is used. 

1. ) For V = we let V = 0. 

2. ) Let V = UW and U ^ ^ W. We have 

pnf (EWiu) = pnf(U a ) pnf (aWw) 
= pnf (U ff ) pnf (Waw) 
= pnf(C/ CT ) pnf (Wr) pnf((a + t)w) 

for some a £ R(U) and r G R(W). Thus, for each p we have: 

V p = min { U a W T | p = a + r, a G R{U),t G R{W) } . 

The minimum search is a little tricky, because U^s may have 
different length so that we are forced to scan through the word 
W T . The overall time we need for these comparisons can be 
bounded however by C(s(C/)||W||) by similar techniques as used 
above in other proofs. (For polynomial time results such a tuning 
is not necessary, and can be omitted.) For creating smaller tables 
we used inductively the time 0(s(U)(\\U\\ + s(U)) + s(W)(\\W\\ + 
s(W)). We need another term (D(s(U)s(W)) for computing the 
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length |pnf(V r p )| and, in any case, we consider all a E R(U) and 
all r E R(V). 

The time we need to do this for all p £ R(V) is bounded by 
0{s(U)\\W\\+s{U)(\\U\\+s(U))+s(W)(\\W\\+s(W))+s(U)s(W)). 
This is within our time bound, because 

s(U)\\W\\ + s(U)(\\U\\ + s(U)) + s(W0(||W|| + s(W)) + s(U)s(W) 
< (s(U) + s(W))(\\U\\ + \\W\\ + s(U) + s(W)) 
= s(V)(\\V\\ + s(V)). 

3.) Let V = aTUpt. We have 

pnf (aTUptw) = a'T pnf (U(eq + 

= olTU a pnf ((<r + eg + 

for some a E R(U) and some a' with a = ep + a' and |a'| < p 

and |e| < 1. 

Now, 

pnf((er + eq + = 0tpni{pw) 

for some /?' and p such that a + eq + (5 = pq + and < q 
and /9 = /up. 

Thus, for each p we may define: 

V p = min { a'TU a 0t } 

where the minimum is taken over all a E R(U), \a'\ < p, and 
|/3' | < q satisfying: 

a = ep + a and a + eg + (3 = pq + and p = //p. 

The time we need to do this for all p E R(V) is bounded by 

0(l + s(U)(\\U\\+s(U)) + s(V)). 

This is within our time bound, because s(U) = s(V) and ||V|| > 
2 + \\U\\. 

□ 

Corollary 7.6 Lei p \ q. We can design a quadratic time algorithm 
running in time C(s(V)(||V|| + s(V)) which solves the following prob- 
lem. 

Input: A valley v 

Problem: Compute the peak normal form pnf(-u). 
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Proof. The Britton reduction of a valley is a valley and can 
be computed in linear time, keeping the integers in binary notation. 
Again in linear time we find a standard valley V and 7 € Z such that 
v ~ 1/7. Theorem 17.51 yields 

pnf (v) = min { V p llnf (p + 7) | p G fl(V) } . 

By Proposition 14.41 llnf (p 4- 7) can be computed in linear time. 

The final step takes time 0{\\V\\ s(V)) because \R(V)\ £ 0{s{V)). 

□ 

7.4 Proof of Theorem 17.11 

We are now in a position to complete the proof of the main Theo- 
rem (7TTJ Let w € {a, A,t,T }* be the input word. First we rewrite w 
as a word in Z * { t, T }* and compute its Britton reduction u in linear 
time. We have < and s(u) < s(w). Next, u is partitioned into 
a\t ■ ■ ■ otktDTfli ■ ■ ■ Tf3 m , where D is horocyclic or difficult. Again, we 
have \\D\\ < \\u\\ and s(D) < s(u). In any case, the peak of u is inside 
D, so, using the reduction given in Theorem 16.14 it is sufficient to 
compute the peak normal form of a constant number of horocyclic or 
difficult words D' with \\D'\\ < \\D\\ + r and s(D') = s(D) and r is a 
constant. From that we get the peak normal form of u in linear time. 

If D' is horocyclic, the peak normal form equals the length-lexico- 
graphical normal form which can be found in linear time, as demon- 
strated in Proposition 14.41 

For the case where D' is difficult, we have shown in Section 17.11 
that we can reduce the computation of the peak normal form to the 
computation of two valleys. The reduction takes linear time. Finally 
peak normal forms for all valleys under consideration can be found in 
quadratic time 0(s(w)(||w|| + s(w)) using Corollary 17.61 □ 

Conclusion 

We have seen how to compute the geodesic peak normal form in 
quadratic time in the case p \ q. Actually, we have shown that the uni- 
form problem is decidable in polynomial time. The uniform problem 
takes as input a word w € { t, t -1 , a, a" 1 }* and two integers p and q 
written in unary with p \ q. As the values r defined in Eqs. [2] and [5] 
are polynomial in q, our quadratic time non-unifom algorithm yields 
a polynomial time uniform algorithm. 

For the general case, where p does not divide q, it remains to com- 
pute geodesies for difficult words, i.e., those, for which the t-sequence 
of the Britton reduction starts with t~ l and ends with t. 
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It is open whether the remaining problem is in P or whether it is 
Co-NP complete or whether the complexity is somewhere in between. 
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Appendix 



Example: Computing length-lexicographical nor- 
mal forms of horocyclic elements 

We give an example to illustrate how our algorithm for finding the 
length-lexicographical normal form of horocyclic elements works. Sup- 
pose we have p = 1, q = 3 and we want to find the length- lexicographi- 
cal normal form of the word (of length 63) 

w = aW W W^a W 3 
= a 7 ta u TA 2 t 2 a 9 Ta 2 Ta 23 
= 7tUT-2tt9T2T23. 

We proceed in several steps. 

Precomputation and finding a suitable slope 

Three Britton reduction steps yield the integer w ~ 157. The greedy 
algorithm of Proposition[0]gives w ~ t52Tl ~ t 2 17TlTl ~ t 3 5T2TlTl, 
hence Unf(u;) = i 3 llnf(u), for the slope u = 5T2T1T1. 

Computing the length-lexicographical normal form for a 
slope 

We now regard the slope u = 5T2T1T1. Since p = 1 and q = 3, we 
have r = 4. We use the algorithm described in the proof of Proposi- 
tion 14.41 Table 17.41 shows the intermediate results of the unoptimized 
algorithm as presented in the first half of the proof. Note that the 
algorithm never has to store more than two columns of the table at 
the same time. So, due to the length of the entries, linear space is 
needed. 

Now let's consider the optimized variant from the second part fo 
the proof. The computational steps are essentially the same, but only 
a constant amount of data is stored. For example for % = 1 we have 
\p\,p\ = min{llnf (u(l, r)) : t € R} = 5, so after having computed the 
column llnf(u(l, p)) in Table [731 we can cut on? the first 5 letters of 
every entry, save their lexicographical order as an order on R (—4 < 
-3 = -2 = -1 < = 1 < 2 = 3 = 4), and output them. The 
complete output of the algorithm is shown in Table 17.41 

Following the pointers and reading off the result, starting at (p, i) = 
(1,3), gives llnf(u) = t2TTT-lT-2. Thus, we have 

llnf(to) = t 3 llnf(u) = t 4 2TTT-lT-2 = iVt _3 a -1 t -1 al -2 . 
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Table 1: Intermediate resulte during the computation of llnf(5T2TlTl) 
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Table 2: Output of the optimized algorithm for llnf (5T2TiT±) 
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